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Abstract 

The extrapolation of small-cluster exact-diagonalization calculations is 
used to examine ferromagnetism in the one-dimensional Hubbard model with 
long-range and correlated hopping. It is found that the correlated hopping 
term stabilizes the ferromagnetic state for a wide range of electron interac- 
tions U and electron concentrations n. The critical value of the interaction 
strength U c above which the ferromagnetic state becomes stable is calculated 
numerically and the ground-state phase diagram of the model is presented 
for physically the most interesting cases. 



The Hubbard model has become, since its inception pQ in 1963, one of the most 
popular examples of a system of interacting electrons with short-range interactions. 
It has been used in the literature to study a great variety of many-body effects in 
metals, of which ferromagnetism, metal-insulator transitions, charge-density waves 
and superconductivity are the most common examples. Of all these cooperative 
phenomena the problem of ferromagnetism in the Hubbard model has the longest 
history. Although the model was originally introduced to describe the band ferro- 
magnetism of transition metals, it soon turned out that the single-band Hubbard 
model is not the canonical model for ferromagnetism. In fact the existence of 
saturated ferromagnetism has been proven rigorously only for very special lim- 
its. The first well-known example is the Nagaoka ferromagnetism that comes from 
the Hubbard model in the limit of infinite repulsion and one hole in a half-filled 
band |2j. Another example where saturated ferromagnetism has been shown to 
exist is the case of the one- dimensional Hubbard model with nearest and next- 
nearest-neighbor hopping at low electron densities [3]. Moreover, several examples 
of the fully polarized ground state have been found on special lattices (special con- 
duction bands) as are the fcc-type lattices j3] , the bipartite lattices with sublattices 
containing a different number of sites [3], the lattices with unconstrained hopping 
of electrons jB] and the flat bands [7j. This indicates that the lattice structure and 
the kinetic energy of electrons, i.e., the type of hopping play an important role in 
stabilizing the ferromagnetic state. 

In our previous paper [H] we have shown that if the electron hopping is described 
by a more realistic model (than the nearest-neighbor hopping) then ferromagnetism 
comes from the Hubbard model naturally for a wide range of the model parameters. 
No extra interactions terms should be included. In particular, we have found that 
the long-range hopping with power decaying hopping amplitudes ty given by [H| 
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gives rise to ferromagnetism for electron densities above half-filling. As soon as q 
that controls the effective range of the hopping (0 < q < 1) is different from zero 
the ferromagnetic state is stabilized for all Coulomb interactions U greater than 
some critical interaction strength U c (q) that value is dramatically reduced with 
increasing q. From this point of view one of main reasons why the ferromagnetic 
state absent in the ordinary Hubbard model with nearest-neighbor hopping (q = 0) 
is that the description of the electron hopping was too simplified. 

Here we further generalize this model by introducing the correlated hopping 
term, in which the cx-electron hopping amplitudes between lattice sites i and j 
depend explicitly on rii_ a and rij- a occupancy, i.e. 



The importance of the correlated hopping term on the ground-state proper- 
ties of the Hubbard model has been already mentioned by Hubbard p. Later 
Hirsch ^0] pointed out that this term may be relevant in explanation of super- 
conducting properties of strongly correlated electrons. Here we examine effects of 
this term on the stability of the fully polarized ferromagnetic state. The same sub- 
ject has been studied recently by Amadon and Hirch as well as by Kollar and 
Vollhardt , however they considered hopping only between the nearest-neighbor 
sites. 

The selection of hopping matrix elements in the form given by Eq. 1 and Eq. 2 
has several advantages. It represents a much more realistic type of electron hopping 
on a lattice (in comparison to nearest-neighbor hopping), and it allows us to change 
continuously the type of hopping and thus immediately study the effect of long- 
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range and correlated hopping on the ground state properties of the Hubbard model. 

The Hamiltonian of the single-band Hubbard model in which the effects of 
long-range and correlated hopping are incorporated is given by 

H = 12 ta i3 c tcja + UY, n^n ib (3) 

ija i 

where cf a and Cj CT are the creation and annihilation operators for an electron of spin 
a at site i, n ia is the corresponding number operator (N a = J2i n ia) an d U is the 
on-site Coulomb interaction constant. 

To examine possibilities for existence of ferromagnetism in this model the 
ground states are determined by exact diagonalizations for a wide range of model 
parameters (q,t',U,N = J2aN a )- Typical examples are then chosen from a large 
number of available results to represent the most interesting cases. The results ob- 
tained are presented in the form of phase diagrams in the U-q plane. To determine 
the phase diagram in the U-q plane (corresponding to some i! and N) the ground 
state energy of the model is calculated point by point as functions of q and U. Of 
course, such a procedure demands in practice a considerable amount of CPU time, 
which imposes severe restrictions on the size of clusters (L) that can be studied 
with this method [L ~ 16). Fortunately, we have found that the ground-state 
energy of the model depends on L only very weakly (for a wide range of the model 
parameters) and thus already such small clusters can describe satisfactorily the 
ground state properties of the model. 

The results of our small-cluster exact-diagonalization calculations obtained on 
finite clusters up to L = 16 sites are summarized in Fig. 1 and Fig. 2. There is 
shown the critical interaction strength U c , above which the ground state is ferro- 
magnetic, as a function of q for selected values of n = N/ L and t' (n = 3/2, 7/4; t' = 
0,0.2,0.4). To reveal the finite-size effects on the stability of ferromagnetic do- 
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mains, the behavior of the critical interaction strength U c (q) has been calculated 
on several finite clusters at each electron filling. It is seen that finite-size effects on 
U c are small and thus these results can be satisfactorily extrapolated to the ther- 
modynamic limit L — > oo. Our results clearly demonstrate that the ferromagnetic 
state is strongly influenced by correlated hopping (£') and generally it is stabilized 
with increasing t'. The effect is especially strong for intermediate and strong values 
of q. Even, there exists some critical value of q above which the ground state is 
ferromagnetic for all nonzero U. With increasing t' this critical value shifts to lower 
values of q (that represent a much more realistic type of electron hopping) and the 
ferromagnetic domain correspondingly increases. Performing exhaustive numerical 
studies of the model for a wide range of electron concentrations (on different lattice 
clusters) we have found that the model exhibits the same behavior for all electron 
concentrations above half- filling [T^], and that with increasing concentration this 
effect becomes more pronounced (see Fig. 1 and Fig. 2). These results clearly show 
that ferromagnetism comes naturally from the Hubbard model with long-range 
and correlated hopping for a wide range of model parameters without any other 
assumptions. This opens a new route towards the understanding of ferromagnetism 
in the Hubbard model. 

In summary, the extrapolation of small-cluster exact-diagonalization calcula- 
tions was used to examine ferromagnetism in the one-dimensional Hubbard model 
with long-range and correlated hopping. It was found that the correlated hopping 
term stabilizes the ferromagnetic state for a wide range of electron interactions U 
and electron concentrations n. The critical value of the interaction strength U c 
above which the ferromagnetic state becomes stable was calculated numerically 
and the ground-state phase diagram of the model is presented for physically the 
most interesting cases. 
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Figure Captions 



Fig. 1. The critical interaction strength U c as a function of q calculated for 
different t' and L at n = 3/2. Curves from up to down correspond to: t' = 0,0.2 
and 0.4. 

Fig. 2. The critical interaction strength U c as a function of q calculated for 
different if and L at n = 7/4. Curves from up to down correspond to: t' = 0,0.2 
and 0.4. 
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